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Abstract 

A detailed derivation of the collisional widths of collective vibrations is 
presented in both quantal and semi-classical frameworks by considering the 
linearized limits of the extended TDHF and the BUU model with a non- 
Markovian binary collision term. Damping widths of giant dipole and giant 
quadrupole excitations are calculated by employing an effective Skyrme force, 
and the results are compared with GDR measurements in Lead and Tin nuclei 
at finite temperature. 
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I. INTRODUCTION 



Excitation energy dependence of the giant resonance is still one of the open problems in 
the study of collective modes in nuclei at finite temperatures. Recent experimental investiga- 
tions of the giant dipole resonance in the mass region ^4=110-140 and the ^osp^ nucleus 
1^ show that damping monotonically increases with excitation energy, and in the former 
case, saturates at high excitations. In medium-weight and heavy nuclei at relatively low 
temperatures the owerwhelming contribution to damping arises from the spreading width 

due to mixing of collective states with more complicated states, which is dominated by 
the coupling with 2p-2h doorway excitations There are essentially two different ap- 

proaches for calculation of the spreading width F^: (i) Coherent mechanism due to coupling 
with low-lying surface modes which provides an important mechanism for damping of giant 
resonance in particular at low temperatures (ii) Damping due to the coupling with 

incoherent 2p-2h states which is usually referred to as the collisional damping HJlO[| and the 
Landau damping modified by two-body collisions |[TT|JT^ . Calculations carried out on the 
basis of these approaches are partially successful in explaining the broadening of the giant 
dipole resonance with increasing temperature, but the saturation is still an open problem 
p!3| . In this work, we do not consider the coherent contribution to the spreading width 
due to the coupling with low-lying surface modes, but investigate in detail the collisional 
damping at finite temperature due to decay of the collective state into incoherent 2p-2h 
excitations in the basis of a non-Markovian transport approach. 

Semi-classical transport models of the Boltzmann-Uehling-Uhlenbeck (BUU) type are 
often employed for studying nuclear collective vibrations [Q. Although, these models give 
a good description for the average resonance energies, they do not provide a realistic de- 
scription for the collisional relaxation rates. In these standard models, the collision term is 
treated in a Markovian approximation by assuming that the two-body collisions are local in 
both space and time, in accordance with Boltzmann's original treatment. This simplifica- 
tion is usually justified by the fact that the duration of two-body collisions is short on the 
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time scale characteristic of macroscopic evolution of the system. However, when the system 
possesses fast collective modes with characteristic energies that are not small in comparison 
with temperature, the standard Markovian treatment is inadequate. It leads to an incorrect 
energy conservation factor in the collision term, which severely restricts the available 2p-2h 
phase space for damping of the collective modes. Therefore, for a proper description of the 
collisional relaxation rates, it is necessary to improve the transport model by including the 
memory effect due to finite duration of two-body collisions [[15| -[T7[|. 



Recently, we have investigated the collisional damping of giant resonances within the 
linearized limit of the BUU model with a non-Markovian collision term, and derived closed 



form expressions for damping width of isoscalar and isovector collective vibrations ||T8|,|T9 
Also, the model has been applied to study the density fluctuations and the growth of insta- 
bilities in the nuclear matter within the stochastic Boltzmann-Langevin approach pO| , pl 



As a result of the non-Markovian structure of two-body collisions, in expressions of transport 
coefficients of collective modes (i.e. damping width and diffusion coefficient), the available 
phase space for decay into 2p-2h states is properly taken into account with the correct energy 
conserving factor. In nuclear matter limit and for isotropic nucleon-nucleon cross-sections, 
by employing the standard approximation familiar in Fermi Liquid theory, it is possible to 
give analytical expression for the collisional width as F = Fo[(/ic<j)^ + (27rT)^], where Fq is dif- 
ferent for different resonance and determined by nuclear matter properties and in-medium 
cross-sections. The quadratic temperature dependence fits well with the measured GDR 
widths in ^"^^Sn and '^^'^Ph nuclei, however the factor Fq calculated with a cross-section of 
40 mb underestimates the data by a factor of 2-3. In a recent work, we have performed 
calculations by employing energy and angle dependent free nucleon-nucleon cross-sections 
and by taking surface effects into account |^ . 



In this work, we give a brief description of the non-Markovian extension of the nuclear 
transport theory in both quantal and semi-classical frameworks, and present a detailed 
derivation of the collisional widths of collective vibrations. The derivation is carried out in 
both quantal and semi-classical frameworks by considering linearized limits of the extended 



TDHF and the BUU model with memory effects. The major uncertainty in the calculation 
of coUisional widths arises from the lack of an accurate knowledge of the in-medium cross- 
sections in the vicinity of Fermi energy. The Skyrme force provides a good description 
of the imaginary part of the single particle optical potential and its radial dependence in 
the vicinity of Fermi energy [^. Therefore, it may be suitable for describing in- medium 
effects in the collision term around Fermi energy. By employing an effective Skyrme force, 
we perform calculations of the damping widths of the giant dipole and giant quadruple 
excitations in semi-classical approximation and compare them with the GDR measurements 
in Lead and Tin nuclei at finite temperature. 



II. ONE-BODY TRANSPORT MODEL WITH MEMORY EFFECTS 

In the extended TDHF approximation, the evolution of the single particle density matrix 
p(t) is determined by a transport equation p^ , 

zh^^p-[h,p]=Kip) (1) 

where h{p) is the mean-field Hamiltonian and the quantity on the right-hand-side represents 
a quantal collision term, which is specified by the correlated part of the two-particle density 
matrix as K{pi) = tr2[v,Ci2] with v as the effective residual interactions. The correlated 
part of the two-particle density matrix is defined as Cu = P12 — P1P2, where piP2 represents 
the antisymmetrized product of the single-particle density matrices, and it is determied by 
the second equation of the BBGKY hierarchy. Retaining only the lowest order terms in the 
residual interactions, the hierarchy can be truncated on the second level, and the hence the 
correlated part of the two-particle density matrix evolves according to, 

ih^C,2-[h,C,2]=Fu (2) 

where 

F12 = (1 - Pl)(l - P2)VP1P2 - P1P2V{1 - Pl)(l - p2)- (3) 
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An expression for the collision term can be obtained by formally solving Eq.(2), and substi- 
tuting the result into Eq.(l), 

K{p) = dTtr2[v,G{t,t - T)Fi2{t - T)G\t,t - t)] (4) 



where 



G{t, t-T)=T- exp[-^ f dt'hit')] (5) 

denotes the mean-field propagator. 

The transport equation (1) is usually considered in semi-classical approximation. In this 
case, one deals with the phase-space density /(r, p) defined as the Wigner transform of the 
density matrix, 

/('''P) = / TTT^e " <P+ oIpIp-TT > • (6) 



Performing the Wigner transform of Eq.(l) and retaining the lowest order terms in gradients 
in accordance with the standard treatment |2^J2^ , yields a semi-classical transport equation 



for the phase-space density, 

|/(r,p)- W),/(r,p)} = W). (7) 

Here, the left hand side describes the Vlasov propagation in terms of the self-consistent 
one-body Hamiltonian h{f), and K{f) represents the collision term in semi-classical ap- 
proximation, which has a non-Markovian form due to the memory effects arising from the 
finite duration of two-body collisions, 

KU) = J dp2dpsdp, I* dTw{12-34;T)[{l - f){l - f2)f3U - fUl - f3){l - U)]t-r (8) 

In this expression the phase-space density is evaluated at time t — r according to fj (t — r) = 
/(r — rpj/m, pj + rVf/, t — r) and the collision kernel is given by 

w{12-34;t) = ^-1^(12; 34)[^i(r)^2(r)(?^(r)(?:(r) + c.c] (9) 



where gjir) is the Wigner transform of the mean- field propagator 

and W{12] 34) denotes the basic two-body transition rate 

^^^2' = (2^1 < ^^l-l^^ >A + p. - P3 - P4) (11) 

which can be expressed in terms of the scattering cross-section as 

The colhsion term involves different characteristic time scales including the average du- 
ration time of two-body collisions t^, the characteristic time associated with the mean-field 
fluctuations Tmj and the mean-free-time between two-body collisions r^. The range of the 
integration over the past history in the collision term (8) (and also in (4)) is essentially deter- 
mined by the duration time of two-body collisions. Usually, two-body collisions are treated 
in a Markovian approximation by assuming the duration time of collisions is much shorter 
than the other time scales <^ Tm/,TA, which would be appropriate if two-body collisions 
can be considered instantaneous. In this case, the r dependence of the phase-space density 
in the collision term can be neglected fj{t — r) ^ fj{t), and the mean- field propagator can 
be approximated by the kinetic energy alone, gj ~ expl—irej] with ej = p|/2m. This yields 
energy conserving two-body collisions, and the resultant semi-classical transport equation is 
known as the BUU model ||2^ . The standard description provides a good approximation at 
intermediate energies when the system does not involve fast collective modes, since the weak- 
coupling condition is well satisfied due to relatively long mean-free-path of nucleons. When 
the system possesses fast collective modes, for example high frequency collective vibrations 
or rapidly growing unstable modes, the Markovian approximation breaks down and the in- 
fluence of the mean-field fluctuations in the collision term becomes important. The finite 
duration time allows for a direct coupling between two-body collisions and the mean-field 
fluctuations, which strongly modifies the collisional relaxation properties of the collective 
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modes as compared to the Markovian limit, in particular at low temperatures |2^. In this 
work we consider the mean-field dominated regime in which the nucleon mean-free-time is 
long as compared to the characteristic time associated with the mean-field fiuctuations and 
the duration time of collisions r^, r^/ ^ ta, which maybe referred as weakly non-Markovian 
regime. In this case, the r dependence of the phase-space density in the collision term can 
be neglected, as before, and the collision term takes essentially a Markovian form with an 
effective transition rate given by Jq drw (12; 34:] t). When all different time scales are of the 
same order of magnitude, the collision term becomes strongly non-Markovian, and the time 
evolution of the system is accompanied by off-shell two-body collisions . 



III. TRANSPORT DESCRIPTION OF COLLECTIVE VIBRATIONS 

We apply the non-Markovian transport model developed in the previous section to de- 
scribe small amplitude collective vibrations around a stable equilibrium in the linear re- 
sponse approximation, and present an explicit derivation of the expression for the collisional 
damping widths of the collective modes in both quantal and semi-classical frameworks. 



A. Quantal treatment 

To describe small amplitude collective vibrations around a finite temperature equilibrium 
state po, we linearize Eqs.(l) and (2) for small deviations 6p = p — po and 6C12 = C12 — C°2, 

d 

i/i— 5p - [5h, Po] - [ho, 5p] = tr2[v, 6C12] (13) 

and 

ih^^SCu - [Sh, C?2] - [ho, SCu] = SF,2 (14) 

where 6h = {dU / dp) q- 6 p represent the small deviations in the single-particle density matrix, 
the quantity 6F12 is 
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SF^2 = -Sp^{l - pDvpIp^, - (1 - p'DSp.vpIp^, + 
(1 - p?)(l - pI)v6^p', + (1 - p?)(l - p°)t;p?5p2 - h.c. (15) 

and the equilibrium correlation function C^^ is determined by 

-[ho,C',,]=F^, (16) 

with F°2 &s the equilibrium value of F^. 

We can analyze the collective vibration by expanding the small deviations 6p in terms 
of normal modes of the system p9[ , 



6p{t)=J2[zx{t)p\ + mp,] (17) 

where p\ and px represent the normal modes of the system. When the damping width is small 
as compared to the mean frequency of the mode, we can follow a perturbation approach and 
determine the normal modes by the standard RPA equations without including the collision 
term, 

nuJxp\-[hlpo]~[ho,pi] = 0. (18) 

Here, uj\ is the frequency of the normal mode and h\ represents the positive frequency part 
of the vibrating mean-field. It is convenient to introduce the RPA amplitudes 0\ and 0\ 
associated with normal modes according to p\ = [0\, po] and its hermitian conjugate. In 
the representation which diagonalizes po, the RPA amplitudes can be expressed as 



, At, < n\h\\m > 

<nO\m>=- ^ 19) 

hujx - e„ + 



and they are normalized as tr[0\,0\]po = 1. Substituting the expansion (17) into Eq.(13) 
and projecting by Ox yields 

^ + iuJxZx = "^^A^A (20) 

for the amplitudes of the normal modes. These amplitudes execute a damped harmonic 
motion with a damping coefficient given by 
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T^ = tr[d^,v]Ci (21) 

and it describes the spreading width of the RPA mode due to couphng with the two particle- 
two hole states, which is usually referred to as the collisional damping. In this expression, 
C\ denotes the positive frequency part of the correlations and it is determined by 

hu^^Ci-[hlcU-[h,,Ci]=Fl (22) 

where F\ represents the positive frequency part of 5Fi2. In the representation diagonalizing 
Po, the correlation can be expressed as 

< nm\[h{,C^2] + Fl\kl > 



< nm\Cl\kl >-- 



huJx -en-em + ek + ei-ii] 
< nm\[dl F^^] + Fl\kl > 



(23) 



hux - en - em + Gk + ei - ir] 
According to expression (16), the matrix elements of the equilibrium correlation is given by 
< nm\Ci2\kl >=< nm\F^2\ki > / {^k + — — — i^)^ which only the principal value 
part is non-vanishing. The second line of the above expression is obtained by replacing the 
energy factors in the intermediate states according to 5{huj\ — e„ — e.^ + eu + ei) and using 
the definition of the RPA amplitudes. Furthermore following the observation 

Pi = -Fl + (1 - p?)(l - pDPI - ^PIPI v]{l - p?)(l - p°) (24) 

the correlation can be expressed as, 

< nm\Cl\kl >=< nm\[dlv]\kl >a IPuPMi - PkRiPuPm] ^35) 
' ^' ^ ^' hux-en-em + ek + ei-i7] ^ ^ 

where p„ denotes the Fermi-Dirac occupation factor, pn = 1 — Pn and < r;,m| [0|, f] >a 
represents the anti-symmetric matrix elements. As a result, the damping width of the RPA 
states is given by |1T6| , 



= 2 H I < nm\[Ox,v]\kl >a \'^5{huJx - e„ - e„ -f- + e,) 

[PkPlpnPm - PnPmpkpl] ■ (26) 
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The same expression for the damping width has been derived in ref. [^] by employing a 
different approach. It also can be obtained in using the Green's function method [^], or 
more intuitive approaches |^. It should be mentioned that the expression (26) was written 
down for the first time by Landau, and it has become a classical result of Fermi Liquid 
theory. However in the nuclear physcis literature several expressions of (26) exist which are 
at variance with Landau's result [^]. The subtlety hinges on the expression (19) for d{. 
Indeed, contrary to the ordinary RPA amplitudes, which at zero temperature have only ph 
or hp components, no phase space factors appear in (19), thus allowing non-zero values of 
< n|0||m > also for pp and hh configurations (see also expression (30), below). 



B. Semi-classical treatment 

It is possible to describe the collective vibrations in semi-classical approximation. In 
this case, one considers the phase-space density 6f{r, p) associated with small amplitude 
vibrations. The equation of motion of the small amplitude vibrations in the semi-classical 
limit is obtained by linearizing the transport equation (7), 

-Sf{r, p) + V • V5/(r, p) - v V5/i-/o = <5K(r, p) (27) 

where 6K{r, p) denotes the linearized collision term, 

6K{r, p) = --J -^—^e-'^dp2 < P + 3 ' P^l b> |p " 3 ' > (28) 

V = p/m and the equilibrium state /(e) is taking to be homogeneous for simplicity. In a 
manner similar to quantal treatment, the phase-space density can be expanded in terms of 
normal modes as 

5/(r, p) = Y^[-ihz^{t)w ■ vol + thzlv ■ VOa]^/ (29) 

where and Ox are the Wigner transform of the RPA amplitudes 0\ and Ox- In the 
perturbation approach, these amplitudes are given by 
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and its complex conjugate. In a similar manner, we can expand the correlation function in 
terms of the normal modes as 6Ci2{t) = J2[zx{t)C\ + zl{t)C\]. By inserting this expansion 
into Eq.(28), we obtain an expression of the collision term in terms of the RPA amplitudes, 

6K{r, p) = ^ E / dp2dpsdp4W{12; 34)[zxAOl + 2*A0a] (31) 

A 

[6{hujx - Ae) - 6{huj, + Ae)][/i/2/3/4 - /1/2/3/4] 

where, Ae = + e, - e, - €3, AOa = Oa(3) + 0^(4) - 0^(1) - 0^(2) with e(j) = pp2m 
and Ox{j) = OA(r, pj), and 1^(12; 34) is the transition rate given by eqs.(ll) or (12). Sub- 
stituting the normal mode decomposition of the phase-space density into (27) and carrying 
out a projection with Ox, we find the expression 

^ 2 (27r^)3 

[S{hujx - Ae) - 6{hujx + Ae)]hhhh 



j dYdpidp2dpzdp4\AOx\'^W{l2;U) (32) 



for the collisional width, where the normal modes in the semi-classical approximation are 
normalized according to 

We note that, this result for the collisional width can directly be obtained by evaluating 
the quantal expression (26) in the Thomas- Fermi approximation [Q. We also note that in 
order to obtain the expressions (26) and (32) for the collisional width in quantal or semi- 
classical forms, the non-Markovian collision term should be linearized in a consistent manner 
by including the contributions arising from the mean-field propagator and the phase-space 
factors. The result is, then, consistent with the Landau's expression for damping of zero 
sound modes, and also, is in accordance with the quantal fluctuation-dissipation relation 
20| , |21|] . If the term involving the mean-field fluctuations (the second term in the left-hand- 



side of eq.(14)) is ignored, one obtains a wrong expression for the collisional damping which 



gives a value that is factor of three larger than its correct value in the nuclear matter [|1^ . 
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It is more convenient to express the semi-classical RPA modes in terms of real functions 
Qx{r,p) and Px{r,p) defined as 

= T^T^^^^""' ^ ^^^^ 

and 

P,(r, p) = . Y^[0*(r, p) - 0,(r, p)]. (35) 
As a result, the normal mode expansion (29) becomes, 

'^/(r,p) = j:M)xl{r,p) +PA(t)x^A(r,p)](-|/) (36) 

where xi — • VPa and Xa — • VQa represent the distortion factors of the phase- 
space density associated with the real variables qx = :y^[z*x + zx] and px = i\f^[zx ~ ^a], 
respectively. In the collision term (31), the factor zxAOx + ZxAOx is replaced by 

^aAO* + zlAOx = ^— (v • V)AO* - z^(v • V)AOa (37) 

^A <^A 



hwx 

where the first fine follows from an identity satisfied by the semi-classical RPA amphtudes, 
0^(r, p) = [hx{Tc)+iTi'v ■'WO\{Yy\/fujJx- In order to deduce the equations for the real variables 
qx{t) and Px{t), we substitute the expansion (36) into Eq.(27) and project the resultant 
equation by Ox and Pa, or equivalently by xi a-iid Xa- This gives two coupled equations 
for qx{t) and px{t), which can be combined to yield an equation in the form of a damped 
harmonic oscillator, 

2 , /rAx2\ Ta . 



Qx+{u^i + {^r)Qx--fQx (38) 



where the coUisional width is given by 



Fa = J drdp,dp2dpsdp4{Axlf + {Axlf]W{12; 34) (39) 



(27rh) 

5{hcux - Ae) - S{hujx + Ae) . - 

TZ /1/2/3/4 

Aricux 
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with the distortion factors normahzed according to 

This expression, which is equivalent to the one given by Eq.(32), provides a useful formula 
to calculate coUisional damping in terms of distortion factors of the momentum distribution 
associated with the collective modes. The distortion factors may be determined from the 
RPA treatment, or can be directly parametrized on physical grounds. In practice, only one 
of the factors, xi or which is associated with a distortion of the momentum distribution, 
contributes the coUisional damping. 

Spin-isospin effects in collective vibration can be easily incorporated in the semi-classical 
RPA treatment by considering proton and neutron degrees of freedom separately. The small 
deviations of the phase-space densities 5/p(r, p), 5/n(r, p) of protons and neutrons are deter- 
mined by two coupled equations analogous to Eq.(27). The collision terms in these equations 
involve binary collisions between proton-proton, neutron-neutron and proton-neutron, and 
a summation over the spins of the colliding particles. Observing that in isoscalar/isovector 
modes protons and neutrons vibrate in-phase/out-of phase, 5/p(r,p) = =l=5/„(r, p), we can 
deduce equation of motions for describing isoscalar/isovector vibrations by adding and sub- 
tracting the corresponding equations for protons and neutrons. Carrying out the semi- 
classical RPA treatment presented above, we obtain P;^ = / drTx{r) with 

n(r) = ^ y dp^dp2dpsdp^[Wpp + Wnn + 2W^] (^) ' Zhhhh (41) 

and 

n = ^ / dp^dp2dp^dp4_{W„ + Wnn) (^) + 2V^p„ ]Zhf2hh (42) 

for the coUisional widths of isoscalar and isovector modes, respectively. Here, Nx = 
! drdp{xxf{-§-J) is a normalization, Axa = Xa(1) + Xa(2) - Xa(3) - Xa(4), Axa = 
Xa(1) - Xa(2) - Xa(3) Xa(4), and Z = [5{nujx - Ae) - 5{nujx + Ae)]/;i(x;A. In these expres- 
sions, transition rates associated with proton-proton, neutron-neutron and proton-neutron 
collisions are given by eq.(12) with the corresponding cross-sections 
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and 

f^l =-^X--U^S + l)\<^^;S,T\v\^^;S,T>^\^ (44) 

Vdr^/p^ (27r;i)3 4;i 8^5^^ ^' 2 ' ' ' ' 2 ' ' ' ^ ^ 

where < pi — P2/2; S", r|f |p3 — P4/2; S, T >a represents the fully antisymmetric two body 
matrix element of the residual interaction between states with total spin and isospin S and 
T. The residual interactions v should be understood as an effective density dependent force. 
It can indeed be shown that a reasonable approximation for v is the so-called G-matrix 
29| . Microscopic G-matrices are not very practical for explicit use, and thus we adopt 



below a more phenomenological point of view replacing the G-matrix by one of the more 
recent Skyrme forces. In doing so, we, however, should be carefull, since in the vicinity of 
nuclear surface Skyrme-type forces usually do not match at all to the free nucleon-nucleon 
cross-sections. 



IV. DAMPING OF GD AND GQ EXCITATIONS 

We apply the formulas (41) and (42) to calculate the collisional widths of the giant 
quadruple and dipole modes by parametrizing the distortion factors of the momentum dis- 
tribution in terms of Legendre functions as xq = p'^P2{cos9) and xd = pPi{cos9). In our 
calculations, we employ an effective Skyrme force, which is parameterized as 

V = to(l + xoP.)(5(ri - r2) + |[(5(ri - r^)^^ + kf^5{r^ - v^)] + 

hkf ■ (5(ri - r2)k + |p°5(ri - r^) (45) 
o 

where k = (pi — p2)/2^ represents the relative momentum operator with k is acting to 
right and k/ is acting to left. In the case of the quadruple mode the collisional width 
is determined by the spin-isospin averaged nucleon-nucleon cross-section, {da/dQ)o = 
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[{da / dVt)pp + {da/dQ)nn + '2{da/dfl)pn]/4:. In the case of the dipole mode the only con- 
tribution comes from the spin averaged proton-neutron cross-section, {da/dfl)pn. In terms 
of the effective Skyrme force these cross-sections are given by 



3- - ^ '^1^,2 , 1 ,2\ , ^3 „i2 



4[to(l + xo) + ^(k^ + k/^) + |p"]^) (46) 



and 



||t„(l + x„) + |(k^ + k/=) + l^oj^) (47) 

where k = (pi — p2)/2^ and k/ = (p3 — p4)/2h are the relative momenta before and after 
the binary colision, and m* denotes the effective mass 

1 1 , 2m 1 , N / M / ^ 

[^ + T2T^i^t, + 5t2)p{r)]. (48) 



m*{r) m % IQ 

In the bulk of the nucleus the Pauli blocking is very effective, and hence, the overwhelming 
contributions to momentum integrals in expressions (41) and (42) arise in the vicinity of the 
Fermi surface. We can approximately perform these integrals by employing the standard 
coordinate transformation 



/ (ipic/p2(ip3(ip4(5(pi + P2 - P3 - P4) ■ ■ ■ ~ / PF^:-deide2de3de4 tttt;^^^ ' ' ' (^9) 

J J 2 cosd / 2 

Furthermore, for temperatures small compared to the Fermi energy, T <^ ep-, the energy 
integrals can be calculated analytically using the formula [p4| , p5| , 

/ de,de,de,de,5{hu ± Ae)f\hM, ^ (M^ + (2^T)^ (50) 

J 6 1 — exp(— n,a;/T) 

Then, the bulk contribution to the coUisional widths of the quadrupole and dipole modes 
can be expressed as. 



r 



m*>p^/Q(r)[(M' + (2^r)1 (51) 



Nq 5 
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and 



h 271 



m*Vi?(r)[(M' + (27rr)']. 



(52) 
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Here, p denotes the particle density, p — (4/(27r^)^)(47r/3)p|^, the normahzations are Nq — 
(47r/5) / drm*pp, Nd = (47r/3) / drm*pp, and the quantities /q. Id are given by 



In these expressions, the angular integrals can be performed analj^ically by noting that 
in the vicinity of Fermi surface the momentum dependent terms in the cross-sections can 
be expressed in terms of the standard variables as k • k/ = —kpsm"^ 9/2cos(f) and k^ = 
k/^ — kpSm^9/2, and cos ^3 = (cos^/2)^ — (sin ^/2)^ cos 0. As already mentioned earlier, 
the cross-sections based on Skyrme forces have a strongly erroneous behaviour at very low 
densities. For this reason, we can not use expressions (51) and (52) far out in the surface. 
It is therefore absolutely necessary to develope effective forces which have the correct free 
cross-section limit. For the time being, we develope an interpolation scheme. In the vicinity 
of nuclear surface, p{r) <^ po, the Pauli blocking is not effective. In this case, it is convenient 
to transform the integration variables in (41) and (42) into the total momenta P = pi -|- p2, 
= P3 + P4 and relative momenta q = (pi — P2)/2, q/ = (pa — P4)/2 before and after the 
collision. Due to the energy conservation, the magnitude of the relative momentum after 
the coUision is restricted according to qf — y/q'^ T m*huj. In the tail region for e^, T" <^ ^o;, 
the expressions (41) and (42) may be estimated by omitting the Pauli blocking factors and 
neglecting the g-dependent terms. This gives. 




(53) 



and 




(54) 




n 27r 



pp%{m*hcuf/^lQ{r) 



(55) 



and 
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rr'(r) ^ ^lppU^*nujy/'lnir). (56) 

We define an effective Pauli blocking factor as the ratio of the damping width with and 
without the Pauh blocking factors in expressions (41) and (42), Fx = rA(r)/r^°-^""'*(r), for 
X = Q or D, and parameterize it in the following form, 

^A(r) = l+(^^)'(F,(0)-1) (57) 

where epiO) is the Fermi energy at the bulk corresponding to the central density po, -Fx(O) is 
the effective factor at the bulk and (3 = huJx/2eF{0). As a function of r, the effective blocking 
factor remains essentially constant and equals to its bulk value until the density reaches about 
1/2 of the central density, and then it smoothly goes to one at the surface of the nucleus. 
This form provides a good approximation for the exact calculations of the effective Pauli 
blocking of 2p — l/i excitations in connection with the collisional damping of single particle 
states as reported in [^. Therefore we expect, it provides a reasonable approximation 
for 2p — 2h excitations, and calculate the collisional widths in an approximate manner by 
smoothly joining the bulk contribution with the surface contribution in accordance with the 
approximate blocking factor (57), 

ra ^ /* (rr W (g||)"°° + rrm - (g||)'l) . /<<.r«M (58) 

and 

We determine the nuclear density in the Thomas-Fermi approximation using a Wood- 
Saxon potential with a depth Vq = —44 MeV, thickness a = 0.67 fm and sharp radius 
Rq = Xm A^l"^ fm. We perform the calculations with the SkM force with parameters a = 
1/6, xo = 0.09, to = -2645Me\//m3, ti = mMeVfm^, t2 = -135Me1//m^ and ts = 
15, 595Me\//m^/2_ por the mass dependence of the resonance energies for spherical medium 
mass and heavy nuclei we use the formulas, hu = 64^-^/3 j^eV for GQR and hu = 80^"^/^ 
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MeV for GDR. Figures 1 and 2 illustrate the relative contribution of the damping widths 
of GDR and GQR for a nucleus with A = 120 as function of r. In these figures and also in 
the other figures, the dashed and dotted lines show the result of the calculations with the 
effective mass m* and the bare mass m, respectively. The sharp rise of r(r) in the vicinity 
of the surface is largely due to the effective mass, which is small in the bulk and approaches 
to its bare value at the surface, and to a lesser extend due to the increase of the Skyrme 
cross-section at low densities. For comparison, the results for constant cross-sections of 
(To = 30 mb and apn = 40 mb are shown in the same figures by solid lines. These constant 
cross-sections correspond to a zero range force with a strength to = —SOOMeVfm^ and 
all other parameters are set equal to zero in (45). Figures 3 and 4 show the atomic mass 
dependence of the GDR and GQR widths and comparison with data, respectively. The 
SkM force with the effective mass underestimates the average trend of GDR for medium 
weight and heavy nuclei by about a factor two. The calculations with the bare mass gives 
a better description of the average trend. In the GQR-case, the discrepancy between the 
calculations and the average trend of data is larger than in the GDR-case. In figures 5 
and 6, the measured GDR widths in and ^°^P6 nuclei are plotted as a function of 

temperature, and compared with the calculations performed with the effective mass and 
the bare mass shown by dashed and dotted lines, respectively. The calculations with the 
effective mass provide a reasonable description of the temperature dependence of data, but 
the magnitude of damping is underestimated in both cases. The calculation with the bare 
mass give larger damping, but the damping widths appear to grow faster than data as a 
function of temperature. 



V. CONCLUSIONS 

In the standard nuclear transport models (mean-field transport models and their stochas- 
tic extensions in semi-classical or quantal form) , the binary collisions are treated in a Marko- 
vian approximation by assuming the duration time of collision is much shorter than the 
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mean-field fluctuations and the mean-free-time between collisions, which would be appropri- 
ate if two-body collisions can be considered instantaneous. As a result, the standard model 
provides a classical description of transport properties of collective motion that is valid at 
low frequency-high temperature limit. When the system possesses fast collective modes, the 
standard description breaks down and it is necessary to incorporate memory effect associ- 
ated with the finite duration of binary collisions. This yields a non-Markovian extension of 
the transport description in which the basic transition rate is modified by involving a direct 
couphng between collective modes and two-body collisions. The extended model leads to 
a description of the transport properties of collective modes that is in accordance with the 
quantal fiuctuation-dissipation relation. In this work, we present a detailed derivation of 
the coUisional widths of isoscalar and isovector collective nuclear vibrations in both quantal 
and semi-classical frameworks by considering the linearized limits of the extended TDHF 
and the BUU model with non-Markovian collision term. The standard treatment with a 
Markovian collision term leads to vanishing coUisional widths at zero temperature. Whereas 
in the non-Markovian treatment the coUisional widths are finite and consistent with the 
Landau's expression for damping of zero sound in Fermi liquids. The numerical result of 
the coUisional damping is rather sensitive to the in-medium nucleon-nucleon cross-sections 
around Fermi energy, for which an accurate information is not available. In the present 
investigations, by employing an effective Skyrme force with SkM parameters, we carry out 
calculations of the damping widths of giant quadrupole and giant dipole excitations in semi- 
classical framework, and compare the results with the GDR measurements in ^"^^Sn and 
208p^ nuclei at finite temperatures. In particular for GDR, the magnitude of the coUisional 
damping with the bare nucleon mass is a sizable fraction of the observed damping widths at 
zero temperature, however the effective mass further reduces the magnitude of damping in 
both cases. Aside from the magnitude, calculations are qualitatively in agreement with the 
broadening of GDR widths as a function of temperature in both and '^^^Pb nuclei. 

One of the main aims of the present investigation was to assess how much of the the total 
width of GRE's is exhausted by decay into the incoherent 2p-2h states. The calculations have 
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been performed within the Thomas- Fermi approximation, which is known from independent 
studies to be very rehable for description of the 2p-2h level densities However, our 

results remain semi-quantitative, since in the Thomas-Fermi framework we need the in- 
medium cross-sections locally down to very low densities, i.e. we need cross-sections which 
interpolate correctly between the free space and the medium. At the moment such cross- 
sections are not available (at least not analytically), and thus we were forced to invent our 
own interpolation scheme, which, although reasonable, is subject to some uncertainties. We 
found that a sizeable fraction of F-'^ is accounted for by the incoherent decay. This is the 
case, for instance, for the GDR and also to a lesser extend for the GQR. In addition, we 
found for the GDR that the percentage of the incoherent decay, depending somewhat on 
the nucleon effective mass, strongly increases with tempearture. This finding is not very 
surprizing, since at temperature T > 3 MeV shell effects are absent and the collectivity of 
the vibrational states is strongly reduced. Therefore to a good approximation, a hot nucleus 
can be regarded as a finite blob of a hot Fermi gas. Inspite of this, in particular at lower 
temperatures, the influence of the low-lying collective sates on the damping is missing in our 
description. Also, the question of the saturation of the GDR width has not been adressed 
in this work. Further studies are needed for a quantitative description of the damping of 
nuclear Giant Resonances. 
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Figure Captions: 

Figure 1: The relative contribution to tlie collisional damping widtli of GDR as a function 
of r for A = 120 at zero temperature. Solid, dashed and dotted lines are calculations 
with a constant cross-section o"p„ = 40 mb, with the SkM force with the effective mass 
and with the SkM force with the bare mass, respectively. 

Figure 2: The relative contribution to the collisional damping width of GQR as a function 
of r for y4 = 120 at zero temperature. Solid, dashed and dotted lines are calculations 
with a constant cross-section o"p„ = 30 mb, with the SkM force with the effective mass 
and with the SkM force with the bare mass, respectively. 

Figure 3: The collisional damping width of GDR as a function of mass number A at zero 
temperature. Solid, dashed, dotted lines are calculations with a constant cross-section 
o'pn = 40 mb, with the SkM force with the effective mass, with the SkM force with the 
bare mass, respectively, and the points show the data. 

Figure 4: The collisional damping width of GQR as a function of mass number A at zero 
temperature. Solid, dashed, dotted lines are calculations with a constant cross-section 
Cpn = 30 mb, with the SkM force with the effective mass, with the SkM force with the 
bare mass, respectively, and the points show the data. 

Figure 5: The collisional damping width of GDR in ^"^^Sn as a function of temperature. 
Dashed, dotted lines and points are calculations with the SkM force with the effective 
mass, with the SkM force with the bare mass and data taken from 0, respectively. 

Figure 6: The collisional damping width of GDR in ^osp^ function of temperature. 
Dashed, dotted lines and points are calculations with the SkM force with the effective 
mass, with the SkM force with the bare mass and data taken from |p, respectively. 
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